Characteristics of bound modes in coupled 
dielectric waveguides containing 
negative index media 



Klaus Halterman, J. Merle Elson, and P. L. Overfelt 

Sensor and Signal Sciences Division, Naval Air Warfare Center, 
China Lake, CA 93555 
Abstract: We investigate the characteristics of guided wave modes in 
planar coupled waveguides. In particular, we calculate the dispersion 
relations for TM modes in which one or both of the guiding layers 
consists of negative index media (NIM)-where the permittivity and 
permeability are both negative. We find that the Poynting vector within 
the NIM waveguide axis can change sign and magnitude, a feature that 
is reflected in the dispersion curves. 
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1 Introduction 

With the continual refinement of laser arrays, and optical fibers, the role of coupled 
waveguide systems in integrated optics has magnified in the past few years. The fabrica- 
tion of arrays which use parallel dielectrics as guiding structures has been demonstrated 
in millimeter-wave dielectric circuits. Furthermore, the coupling of waves in planar di- 
electric sheet waveguides provides a means of reflectionless signal transfer from one 
waveguide to another. For these systems, generally the constituent dielectric materials 
are presumed to have positive values of the electric permittivity, e and magnetic perme- 
ability, /i. Recently, however there has been interest in researching materials which have 
the unusual property of both e and fi being negative, a concept first proposed by Vese- 
lago [1] long ago. These materials support electromagnetic wave propagation in which 
the phase velocity is antiparallcl to the energy flow, or group velocity. Furthermore, 
if both e, and \i are negative in a medium, the refractive index has the extraordinary 
property of also being negative, therefore such media are given the name negative in- 
dex media (NIM), compared with the more familiar positive index media (PIM)-where 
both n and e are positive. As discussed by Veselago [1], a multitude of other anoma- 
lous characteristics of NIM follow from Maxwell's equations, including reversal of the 
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Fig. 1. Coupled waveguide structure considered in this paper. The red and blue 
regions correspond to the waveguide channels with dispersive material parameters. 
The permittivity and permeability in the channels can take positive or negative 
values. The free space parameters take the values eo = Mo = 1- 

Doppler shift, counter-directed Cherenkov radiation cone, and the refocusing of EM 
waves from a point source. It was shown [2] that a NIM slab can focus focus both 
near and far fields, attaining perfect resolution due to the amplification of the evanes- 
cent waves. Subsequently however, the original "perfect lens" model was modified [3] to 
overcome limitations of imperfect materials by introducing a multilayer stack. While the 
study of various NIM structures has intensified during the past few years, the number 
of theoretical works involving waveguides with NIM components is less common. The 
S parameters and electromagnetic properties were investigated [4] numerically for a T 
junction NIM waveguide configuration, and a simple plane two-layered waveguide was 
shown [5] to have a slow-wave factor that tends to infinity at small frequencies. 

On the experimental side, the fabrication of devices constructed with NIM is less 
prolific, as there are no naturally occurring NIM in nature, and thus one must artificially 
engineer a composite system having the desired electromagnetic properties. One such 
medium was constructed [6] in which the effective permeability and effective permittivity 
were shown to be simultaneously less than zero over a finite frequency band. These 
structures that simulated a NIM were comprised of a periodic array of unit cells, where 
each cell contained two split metallic rings located on a dielectric panel along with one 
copper wire. The overall medium is considered homogeneous since the elements of each 
cell varied over length scales much smaller than the wavelengths used in the experiments. 
The experimental situation involving these unconventional materials is thus in the early 
stages of development, and there exists a broad range of potential applications, including 
optical and microwave devices. 

It is the purpose of this paper to investigate the electromagnetic properties of a cou- 
pled planar dielectric waveguide system. In order to study the interaction between NIM 
and PIM configurations, our method will allow for a wide range of material parameters. 
The paper is organized as follows. In Sec. 2, we introduce the geometry and method 
used to calculate the relevant quantities of interest. In Sec. 3, we calculate the dispersion 
relations, electromagnetic fields, and energy flow characteristics for differing waveguide 
structures. Finally in Sec. 4, we summarize our results. 



2 Method 



In this section we give a brief overview of the model used for our coupled waveguide 
system. The schematic of the planar dielectric waveguide is illustrated in Fig. 1. The 
structure consists of two waveguide channels embedded in a medium characterized by 
the free space parameters eo, and /j,q. The upper waveguide layer of width 2d\ has 
permittivity e\ and permeability fii, and the lower guide of width 2d^ has permittivity 
£2 and permeability fi2- The two channels are separated by a central film that is 2d in 
width and is divided in the plane at x = 0. The z axis of the coordinate system coincides 
with the axis of the guide. 

We are interested in guided wave modes, which are electromagnetic waves that prop- 
agate along the waveguide with a given phase and group velocity, intensity distribu- 
tion, and polarization. These modes are relevant since the corresponding fields satisfy 
the wave equation throughout the structure, and the appropriate boundary conditions. 
Furthermore, each mode is characterized by its frequency uo and its wave vector k z . To 
proceed, we assume a sinusoidal time dependence exp(— iujt) for the fields. The dielectric 
media is linear and isotropic in the y and z directions. The translational invariance in 
the y and z direction allows, (for the case of TM modes) one to write the magnetic field 
H as, H = yh[x) exp[i(k z z — cut)], where we have factored out the spatial variation of 
the field in the x direction, and k z is the z component of the wave vector. Employing 
Maxwell's equations, one can immediately write down the reduced wave equation that 
must be solved in each of the layers: 

^-A£ + W (w)c-(w)^]fc(aO=0, 1 = 0,1,2, (1) 

with an analogous expression for the E field. In what follows, we consider material 
parameters in the NIM regions to be dispersive, otherwise the energy density would be 
negative [1]. We take the permittivity and permeability to have the form, 

6i(o;) = l-^, M iM = l-4> ( 2 ) 

where u> e , and u> m are the effective electrical and magnetic plasma frequencies [2], re- 
spectively. To simplify the model, we have neglected damping. Elementary solutions to 
Eq. (1) admit the following wave vectors, 



^)^ 2 _ k 2 i = 0>lj2j (3) 



where the wave vectors are labeled according to the respective values of the permeability 
and permittivity in each region depicted in Fig. 1. It is clear that the longitudinal wave 
vector k z , common to all the layers, serves as a connection to constrain the transverse 
components. 

3 Results 

3.1 NIM Waveguides 

We consider first the case where the permittivity and permeability are both simultane- 
ously negative in each waveguide core, and each of the guiding layers has the same width 
di = c?2 (see Fig. 1). To derive the TM modes, we proceed to seek elementary solutions 
of Eq. (1) that do not radiate away from the system. Upon invoking the continuity of 
E z and H y at the interfaces, we arrive at the following eigenvalue equation or dispersion 
relation for the symmetric and antisymmetric modes: 

± e - 2kod sin(2/Mi) - sin(2/Mi - 2tp) = 0, (4) 




o.8 r i 

0.9 1 1.1 1.2 1.3 1.4 1.5 



C0/C0q 

Fig. 2. Dispersion curves for waveguide structure with material parameters given 
by Eqs. (2). The widths of the middle layer and the two waveguide channels are 
given by d = O.lArj, and d\ = 0.2Aq respectively, where loo = 27re/Arj. The red curve 
corresponds to the symmetric modes, and the blue curve is for the asymmetric 
modes. The central green curves correspond to the dispersion relation for a single 
NIM waveguide of thickness 2di . The inset is a magnification of the same quantity 
localized around one of the bending features. The numbers label points of interest 
as discussed in the text. 

where we have introduced a phase angle ip, 

ijj = arctanf 01 
\kie 

We show in Fig. 2, the dispersion relation for the case where all the dielectric layers are 
made of NIM. We plot k z normalized by a design frequency ljq /c versus the dimensionless 
frequency u)/u)q. In the following, we take io e = 2.34o;o, and to m = 1.98ljq. All curves are 
bound by the light line, k z = (w/c) v /e u 7Io, and the limiting line k z = (ui / c) ^ 'ei(w) jJLi(u>) ■ 
Thus, the "phase-space" available for guided wave modes is contained within the area 
spanned by these two lines, and increases rapidly upon decreasing u>, due to material 
dispersion. The figure illustrates the dispersion diagram with a finite number of bands, 
which arise from the geometrical constraint placed upon the fields from the bounding 
surfaces. Since the fields outside the waveguide channels decay over the length scale k^ 1 , 
regions of the dispersion curves near the light line correspond to a larger decay length, 
thus the fields reside mainly outside the guiding layer. The opposite is true for the to and 
k z dispersion pair which reside in the vicinity of the line k z = (lu/ c) \J t\(^S)\x\{ui). The 
effect of increasing the coupling (decreasing d) between the two waveguide cores has the 
trend of increasing the divide between the two modal bands with different symmetry. 
Although not visible in the plot, the zeroth order TM mode is cutoff in frequency, in 
contrast to the usual PIM case where there is no frequency cutoff. By way of comparison, 
we also show in Fig. 2, the dispersion curve for a single NIM waveguide. It is apparent 
that the symmetric and antisymmetric modes are centered about the single waveguide 
case. One of the most dramatic features of the dispersion relation is the bending back of 
the bands close to the light line. This is a behavior seen in general whenever a waveguide 
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channel consisting of NIM is in proximity to a material with PIM characteristics. The 
inset of Fig. 2 shows a close up of the second modal branch. Some salient features to be 
noted are the change in slope of the dispersion curve, varying from positive to negative as 
one cycles through the points 1 — 3 respectively. We thus find some peculiarities in the 
dispersion diagram that clearly deviate from the well understood coupled waveguide 
system with conventional dielectric material. As will be discussed below, there is an 
abrupt reversal of energy flow in the z-direction between the two different media, and 
this quite naturally is reflected in the allowed guided wave modes. In order to further 
investigate this anomalous behavior of the energy flow in the vicinity of a bend, it is 
necessary to calculate the Poynting vector, which requires the full electromagnetic fields. 

Having obtained the set of k z and w that are solutions to Eq. (4) , the spatially varying 
electromagnetic fields can be straightforwardly calculated from Eq.(l). To demonstrate 
this, we write down the magnetic fields consistent with the dispersion relation in Eq.(4), 

h ^= n e~ k ^- L \ x>L, (6a) 

cos(fcidi — ip) 



cos(fci<ii — tp) 
cos(2/cidi — ip) 
cos(fciG?i — ip) cosh(fc rf) 
cos[fci(x + L) — ip] 



cos(fci<ii — ip) 
cos(ip) 



cosh(fcoa;), \x\ < d, (6c) 
-d < x < —L, (6d) 



COs(fciG?i — Ip) 



k (x+L)^ X <-L, (6e) 



where L = 2<ii + d, and for brevity we show only the symmetric field, noting that the 
antisymmetric fields are calculated similarly. The magnetic field distribution in Eq. (6) 
reflects the fact that the fields are oscillatory within the two waveguide cores, and ex- 
ponentially decaying in the outer regions. For a given set of material parameters, Eqs. 
(6) supplemented with Eq. (4) completely determine (apart from an unimportant over- 
all constant amplitude factor) the electromagnetic properties of the structure. We are 
ultimately interested in the time- averaged power flow density however, as given by the 
real part of the Poynting vector: S = [c/87r](E x H*). Upon calculating the electric field 
distribution from Eq. (6) by taking the appropriate derivatives via Maxwell's equations, 
we find that the transverse component, S x is imaginary, and thus does not contribute 
to the time-average flux of energy. The z component however does contribute, and is 
given by, 

S z (x;k Zl u) = ^-^-{h(x)} 2 , i = 0,1,2. (7) 

In this instance, the direction of S z depends upon the sign of the permittivity of the 
medium, and is discontinuous across the boundaries. On the other hand, the longitudinal 
wave vector k z remains continuous across each layer, as required by translational invari- 
ance. To illustrate the spatial dependence of 5*2, we show in Fig. 3, the time averaged 
Poynting vector throughout the waveguide system. Each of the three curves corresponds 
to an allowed k z and uj pair identified by the three points within the inset of Fig. 2. 
The most visible behavior is the rapid decline in S z in the regions characterized by 
material parameters eo an d ^o, while the two NIM waveguide channels show relatively 
little variance. The blue curve reveals that more total energy is flowing in the central 
and two outer layers, while the parameters corresponding to the green curve shows that 
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Fig. 3. The normalized z-componcnt of the Poynting vector S z as a function of 
dimensionless position. The continual decline in amplitude of S z in the central film 
and the outer two layers from blue to green coincides (in order) with points 1-3 in 
the inset of Fig. 2 respectively. Meanwhile, the two waveguides have relatively little 
shift in amplitude. The vertical lines serve as guides to the eye. 



the net negative S z in the waveguide regions is greater overall. This trend continues 
at other points, so that the sign of the slope of the dispersion relation correlates with 
that of the net energy flow. Since cLo/dk z changes sign, there must be points in the 
dispersion curve that result in S z in the middle and outer layers exactly countering S z 
in the two waveguide channels. This is shown by the intermediate red curve in Fig. 4. 
For this special point on the dispersion curve, dui/dk z = (see point 2, inset of Fig. 2), 
and the net S z is zero. 

The interesting properties of the net energy flow in our coupled waveguide system 
can be investigated further by calculating the time average of the power flow P z across 
all the layers. In order to do this, we spatially integrate Eq.(7) over all of the spatial 
region encompassing the system. For the single layer case, P z can be simply expressed 
as 

d a \- [°° a Qi u ^ c 2 fc z di rcos 2 (fcidi) 1 f sin(2fc]di)i\ , Q , 
P z (k z ,u>)= dxS z (x;k z ,uj) = \ — — + — 1 + \\. (8) 

Figure 4 shows the power across the whole system as a function of modal frequency. We 
see that the single layer result is always between the symmetric and antisymmetric cases. 
The positive P z curves correspond to the lower branch in the dispersion diagram, below 
the point where du>/dk z = (see Fig. 2). Likewise, the negative P z curves correlate 
with those u>, k z pairs above that point. Thus a clear relationship exists between the net 
power and the slope of the dispersion curves: P z is always negative for values of k z and 
u) that satisfy dcu/dk z < 0, and is positive when duj/dk z > 0. For those values of k z and 
u> where P z = 0, we again find consistency with the dispersion curves, where at those 
same points dio/dk z = 0. This result follows from the exact cancellation of negatively 
directed energy flow in the waveguides with the positive contributions from the outer 
PIM layers. 
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Fig. 4. The normalized z component of power, P z , is plotted as a function of di- 
mcnsionlcss frequency lo/loq. The color labeling of the curves coincides with the 
dispersion curves displayed in the inset of Fig. 2. For comparison purposes, the 
green curve given by Eq. (8), is multiplied by a factor of 2. 



3.2 NIM/PIM Waveguides 

Next, we consider the case where the two waveguide channels possess opposite material 
characteristics, i.e., t\ = —e 2l and /ii = —^2- After a lengthy calculation, the dispersion 
relation can be expressed as: 



In Fig. 5 wc present the relationship between k z and uj for the present asymmetric waveg- 
uide configuration. The two differently colored curves shown correlate to the choice of 
sign in Eq. (9). Since there is no longer any symmetry in the x direction, the overall be- 
havior of the dispersion curves is more complex than the previous case, where for a given 
frequency range, each waveguide had identical negative permittivity and permeability 
values. The figure illustrates the unusual change in sign of du)/dk z that occurs when 
traversing the two visible U-shaped branches. To proceed with the calculation of the 
Poynting vector, we again calculate the electromagnetic fields within the five different 
layers. We find, 



y/l - e -4feod S in(2fcidi) ± sin(2-0) = 0. 



(9) 



h(x) 




sin[2(fcidi + ip)} fcn(x _ 2d) 



(10a) 



d < x < L, 



(10b) 



sin[2(fcidi)] \ , \x\ < d, (10c) 




cos[fci(:c + L) + , — d < x < —L, 



(lOd) 



x < -L. 



(lOe) 



From the plot in Fig. 5, we select four points (k z ,ui) on the dispersion curves intersected 
by the oj/uiq = 1.24 line. Inserting these points into Eqs. (10) and then into (7) results 
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Fig. 5. Dispersion relation for the coupled waveguide structure consisting of a NIM 
channel and a PIM channel. The widths of the middle layer and the two guide 
channels are given by d = O.lAo, and d\ = 0.2Ao respectively. The red (blue) curves 
correspond to the plus (minus) sign in Eq. (9). 

in the normalized spatially dependent energy flow that is shown in Fig. 6. The upper 
panel in Fig. 6(a) depicts a relatively flat red line, which is actually three curves merged 
together, along with a dominant energy distribution (black curve) for the parameter 
values uj/loq = 1.24, and k z /(ujo/c) = 2.22. Thus a feature arises in which three out of 
the four wave vector values have negligible contributions to S z in this region. Proceeding 
to the lower panel (b), we exhibit S z for the same parameter values as in (a), except 
we consider now the NIM region. In addition to the energy flow reversing direction, 
all four k z and 10 pairs are now visible and have a relatively complicated asymmetry 
when compared to the case when the waveguide structure is physically symmetric, as 
in Fig. 3, where the guiding layers are both NIM media. It is clear therefore, that there 
can be cases where the energy flow is dominant in either the NIM or PIM slab, or is 
shared between both slabs. 

4 Conclusion 

In this paper we have examined the geometrical and material dispersion properties of 
a coupled planar waveguide structure that contains negative index media. We calcu- 
lated the TM guided wave modes and found a variety of unconventional results. For 
the case where both coupled waveguides have NIM characteristics and are embedded 
in PIM media, the dispersion relation revealed clear deviations from the well under- 
stood case involving positive material parameters. We found that due to the reversal 
of energy flow along the propagation direction, the dispersion curve reflected a striking 
behavior in curvature. This followed ultimately from the correlation between the net 
power in the system and the slope du>/dk z . Furthermore, decreasing the relative dis- 
tance between each wave channel has the effect of further splitting the symmetric and 
antisymmetric dispersion curves. Finally, we presented results for a more complicated 
asymmetric waveguide structure comprised of both NIM and PIM guiding layers. In this 
case, the Poynting vector and dispersion characteristics exhibited nontrivial behavior. 
In particular, we found that there are cases when the energy flow reverses direction 
when going from the NIM to PIM layer, and that the relative distribution of power 




_15 r 1 

-0.5 -0.4 -0.3 -0.2 -0.1 



X/Xq 

Fig. 6. The spatial dependence of the (normalized) Poynting vector S z , in both 
waveguide cores, as a function of x/\q. The top panel (a) depicts S z in the PIM 
guide, and the bottom panel (b) corresponds to the NIM guiding layer. Energy flow 
is always negative in the NIM waveguide, and positive in the PIM waveguide. For all 
four curves, uj/loq = 1.24 and k z /(uJo/c) = 1.27 (blue), 1.59(green), 1.89 (red), and 
2.22 (black). Upon comparing both panels, it is evident that for k z /(LJo/c) = 2.22, 
the PIM waveguide layer bears most of the energy flux from the electromagnetic 
fields. 

varied dramatically between the two, depending on the particular guided wave mode 
under consideration. 
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